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Abstract 



We define signed dual volumes at all dimensions for circumcentric dual meshes. We show that for pairwise Delaunay triangulations 
with mild boundary assumptions these signed dual volumes are positive. This allows the use of such Delaunay meshes for Discrete 
^ ^ Exterior Calculus (DEC) because the discrete Hodge star operator can now be correctly defined for such meshes. This operator 
^—^ is crucial for DEC and is a diagonal matrix with the ratio of primal and dual volumes along the diagonal. A correct definition 
requires that all entries be positive. DEC is a framework for numerically solving differential equations on meshes and for geometry 
processing tasks and has had considerable impact in computer graphics and scientific computing. Our result allows the use of DEC 
^Pwith a much larger class of meshes than was previously considered possible. 
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1. Introduction 

Discrete Exterior Calculus (DEC) is a framework for nu- 
merical solution of partial differential equations on simpli- 
cial meshes and for geometry processing tasks [4, 8]. DEC 
has had considerable impact in computer graphics and sci- 
entific computing. It is related to finite clement exterior 
calculus and differs from it in how inner products are de- 
fined. The main objects in DEC arc p-cochains, which for 
the purpose of this paper may be considered to be a vector 
of real values with one entry for each p-dimensional simplex 
in the mesh. For p-cochains a and b their inner product in 
DEC is a T * p b where * p is a diagonal discrete Hodge star 
operator. This is a diagonal matrix of order equal to the 
number of p-simplices and with entries that are ratios of 
volumes of p-simplices and their (n — p)-dimensional cir- 
cumcentric dual cells. 

For this to define a genuine inner product the entries have 
to be positive. Simply taking absolute values or consider- 
ing all volumes to be unsigned does not lead to correct so- 
lutions of partial differential equations. See Figure 1 to see 
the spectacular failure when unsigned volumes are used in 
solving Poisson's equation in mixed form. That figure also 
shows the importance of the Delaunay property and of our 
boundary assumptions and the success of the signed dual 
volumes for such meshes that we describe in this paper. 

When DEC was invented it was known that completely 
well-centered meshes were sufficient but perhaps not nec- 
essary for defining the Hodge star operator. (A completely 
well-centered mesh is one in which the circumccnters are 
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contained within the corresponding simpliccs at all dimen- 
sions. Examples are acute-angled triangle meshes and tetra- 
hcdral meshes in which each triangle is acute and each 
tetrahedron contains its circumcenter.) For such meshes, 
the volumes of circumcentric dual cells are obviously pos- 
itive. For some years now there has been numerical evi- 
dence that for codimension 1 duality some (but not all) 
pairwise Delaunay meshes yield positive dual volumes if 
volumes are given a sign based on some simple rules. More- 
over, these seem to yield correct numerical solutions for a 
simple partial differential equation. A pairwise Delaunay 
mesh (of dimension n embedded in R w , N > n) is one in 
which each pair of adjacent n-simplices sharing a face of 
dimension n — 1 is Delaunay when embedded in R™. (Imag- 
ine a pair of triangles with a hinge at the shared edge and 
lay the pair flat on a table.) This generalizes the Delaunay 
condition to triangle mesh surfaces embedded in three di- 
mensions and analogous higher dimensional settings. For 
planar triangle meshes and for tetrahedral meshes in three 
dimensions pairwise Delaunay is same as Delaunay. 

We give a sign convention for dual cells and a mild as- 
sumption on boundary simpliccs. With these in hand, for 
pairwise Delaunay meshes it is easy to see that the codi- 
mension 1 dual lengths are positive in the most general 
case (dimension n mesh embedded in M. N ). In addition, we 
prove that such triangle meshes embedded in two or three 
dimensions have positive vertex duals and that the duals 
of vertices and edges of tetrahedral meshes in three dimen- 
sions are positive. This settles the question of positivity for 
all dimensions of duality for simplicial meshes used in phys- 
ical and graphics applications and opens up the possibility 
of using DEC with a much larger class of meshes. 

Note that the results of this paper are relevant for the 



assembly of the duals from elementary duals (defined in 
Section 2) . Thus it is important for algorithms which com- 
pute these dual volumes piece by piece from the elementary 
duals (such as are used in the software PyDEC [2]). It may 
be possible to have alternative formulas which bypass the 
assembly process and give the dual volumes directly. Such 
formulas are hinted at in [9]. Zobel [11] has also hinted at 
the positivity of the duals in some cases but without proofs 
or details. 

There are alternatives to the discrete diagonal Hodge 
star. In finite element exterior calculus, the role of discrete 
diagonal Hodge star is played by the mass matrix corre- 
sponding to polynomial differential forms [1] . This matrix 
is in general not diagonal. Sometimes these are referred to 
as the Galerkin Hodge stars. Our results do not apply to 
that case directly. For finite elements, there are other qual- 
ity measures that may be important. Sec for instance [10]. 

Notation: We will sometimes write the dimension of a 
simplex as a superscript. The notation r -< a means that 
t is a face of a. Circumcenter of a simplex r is denoted c T 
and the circumcentric dual of r as *r. Some figure labels 
like -kp stand for dual of a p-simplex. 

2. Signed Circumcentric Dual Cells 

The dual mesh of a (primal) simplicial mesh is often 
defined using a barycentric subdivision. For every p- 
dimensional simplex, there is a dual (n — p)-dimensional 
cell where n is the dimension of the simplicial complex. 
In DEC circumccntcrs are used instead of barycenters 
because the resulting orthogonality between the primal 
and dual is an integral part of the definition of some of 
the operators of DEC [8]. The circumcentric dual cell of 
a p-dimensional primal simplex r is constructed from a 
set of simplices incident to the circumcenter of r. These 
are called elementary dual simplices with vertices being a 
sequence of circumcenters of primal simplices incident to 
t. The sequence begins with the circumcenter of t, moves 
through circumcenters of higher-dimensional simplices a 1 
and ends with the circumcenter of a top-dimensional sim- 
plex a 11 such that r -<; a p+1 -<•••-< a n . Taking each of the 
possibilities for a % at each dimension i yields the full dual 
cell. 

In the past in the software PyDEC the volume of the dual 
cells has been taken to be the sum of unsigned volumes of 
the elementary dual simplices. Instead we define the vol- 
ume of the dual cells as the sum of signed volumes of its 
elementary dual simplices. Our contribution is in defining 
the sign convention and describing with proofs the class of 
meshes for which the dual volumes and hence Hodge star 
entries are positive. If the primal complex is completely well 
centered every elementary dual has a positive volume. In 



general, the sign of the volume of an elementary dual sim- 
plex is defined as follows. Start from the circumcenter of t. 
Let v p be the vertex such that v p * r is the simplex a p+1 
formed by the vertices of r together with v p . Similarly, for 
p + 1 < i < n — 1, let Vi be the vertex such that Vi * a 1 
is the simplex a t+1 . If the circumcenter of a p+1 is in the 
same half space of <r p+1 as v p relative to r, let s p = +1, 
otherwise, s p = — 1. Likewise, for p+l<i<n — 1, if the 
circumcenter of a t+1 is in the same half space as Vi relative 
to <t' 1 , let Si = +1, otherwise, Si = —1. Then, the sign s of 
the elementary dual simplex is the product of the signs at 
each dimension, that is, s — s p s p+ i ■ ■ ■ s n . 

For illustration of this sign rule, we now consider various 
cases in two and three dimensions. The first example is the 
dual of an edge ab in a triangle abc. From the midpoint of ab 
— its circumcenter - we move to the circumcenter of triangle 
abc. If this move is towards vertex c, then the sign is s = 
s\ = +1, but if it is away from vertex c, as it will be if the 
angle at vertex c is obtuse, then the sign is s = Si = — 1. The 
next example is the dual of vertex a in triangle abc. We will 
consider the simplex formed from the circumcenter of a, the 
circumcenter of ac, and the circumcenter of abc. The move 
from a to the midpoint of ac gives si = +1, since vertex c 
and the midpoint of ac are in the same direction from a. 
The move from the midpoint of ac to the circumcenter of 
abc gives s 2 = +1 if we go towards b and s 2 = —1 if we move 
away from b. The sign of the volume of this contribution to 
the dual of vertex a is s = si s 2 = s 2 . 

For a tetrahedron abed we can expand on the cases for 
triangle abc. For the dual to face abc, we move from the 
circumcenter of abc to the circumcenter of abed. If the cir- 
cumcenter of abed is in the same half space as vertex d rela- 
tive to abc, this move is towards d, the sign is s — si = +1, 
and the signed length (volume) is positive; otherwise, it is 
negative. Of the two contributions to the dual of edge ab, 
we focus on the simplex formed from the circumcenter of 
ab, the circumcenter of abc and the circumcenter of abed. 
The sign si is determined as it was for the dual of edge 
ab in triangle abc. The sign s 2 is +1 if vertex d and the 
circumcenter of abed are in the same half space relative to 
abc. Thus for the dual of edge ab, the sign of the volume is 
s = si s 2 , and both si and s 2 can be either positive or neg- 
ative. As a final example, consider the simplex formed from 
vertex a, the circumcenter of ac, the circumcenter of abc, 
and the circumcenter of abed. This simplex contributes to 
the dual of vertex a. Signs si and s 2 are the same as they 
were for the dual of vertex a in triangle abc. Sign s 3 is — 1 
if triangle abc separates vertex d from the circumcenter of 
tetrahedron abed. The sign of this elementary volume then 
is s = si s 2 s 3 . 

The significance of the sign rule defined above is that it 
orients the elementary dual simplices in a particular way 
with respect to the dual orientation for a completely well- 
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Fig. 1. Solution of Poisson's equation — Am = / in mixed form. In mixed form this equation is the system a = — gradu and diver = /. The 
boundary condition is constant influx on left and outflux on right. The correct solution is linear u which varies only along x-direction and 
a constant horizontal a. The top row shows u and bottom row shows a. The first column shows the correct solution using the results of 
this paper on a Delaunay mesh with correct boundary simplices. The next three columns show various failure modes of alternatives. Second 
column is for unsigned duals using the same mesh as first column. The third column has a single bad (i.e., not one-sided - see Section 4) 
boundary triangle shown shaded in a Delaunay mesh. The fourth column is a non Delaunay mesh. 
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Fig. 2. Examples of sign rule application in 2d. The dot marks the 
circumcenter and green and red are used to denote positive and 
negative volumes respectively. 

centered simplex. Consider two rt-dimensional simplices a 
and <j w which have the same orientation but such that a w 
is well-centered. We are given a bijection between the ver- 
tices of these two simplices such that the resulting simpli- 
cial map is orientation preserving. This vertex map induces 



a bijection between faces of the two simplices and between 
their elementary duals. Let r and t w be two correspond- 
ing p-dimcnsional faces in the two simplices and consider 
their duals *r and *r m . If we consider two correspond- 
ing elementary duals in *r and *t w we can affinely map 
these such that the first vertex (the circumcenter of r or 
T w ) is mapped to the origin and the others are mapped to 
+1 or —1 along a coordinate axis. For the elementary dual 
in -kr w we always choose +1 for all n — p coordinate axes. 
For *rwe choose +1 if the sign along that direction of the 
elementary dual is positive according to the sign rule de- 
scribed above and — 1 otherwise. It is clear (and is easy to 
show using determinants) that the orientation of the corre- 
sponding elementary duals will be same if an even number 
of — 1 directions are used for the elementary dual in * r and 
the orientations will be opposite otherwise. Thus we have 
shown the following result. 

Theorem 1 With a, o~ w , t, andr w as above, the orienta- 
tion of*T is same as that of*T w if an even number of — 1 
signs appear according to sign rule and is opposite otherwise. 
If the orientation of *r is same as *t w we will assign a 
positive volume to *t and otherwise a negative volume. 
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3. Signed Dual of a Delaunay Triangulation 

We first consider the codimension 1 case in the most gen- 
eral setting of a simplicial complex of arbitrary dimension n 
embedded in dimension N > n. After that we consider cases 
other than codimension 1 but in more restricted settings. 
For these latter cases we restrict ourselves to the physically 
most useful cases of triangle meshes embedded in two or 
three dimensions (n — 2 and N = 2 or 3) and tetrahedral 
meshes embedded in three dimensions (n = N = 3). We 
conjecture that these results can be extended to the more 
general setting of arbitrary n and N > n but those cases 
are not as important for physical applications and we leave 
those for future work. For the general codimension 1 case 
we first prove the following basic fact about circumcenter 
ordering for Delaunay pairs. 

Lemma 2 (Circumcenter Order) Let t be an (n — 1)- 

dimensional simplex in R™. Let L and R be points such that 
A = L * t and p = R* r form a non-degenerate Delaunay 
pair of n- dimensional simplices with circumcenters c\ and 
c p , respectively. Then, c\ and c p have the same relative 
ordering with respect to t as L and R. 

PROOF. Consider the collection of (n — l)-dimensional 
spheres containing the vertices of r. Since A and p are a non- 
degenerate Delaunay pair, their circumspheres are empty 
and belong to this collection. It is then easy to see that c\ 
and c p will be in the same order as A and p. See Figure 3. 
For an alternative, more algebraic and detailed proof, see 
Appendix A. (In fact there we show the stronger result that 
the correctness of the circumcenter ordering is equivalent 
to the simplices being a non-degenerate Delaunay pair.) 




Fig. 3. For a Delaunay pair the ordering of the circumcenters is the 
same as that of the top dimensional simplices. See Lemma 2. 

The above lemma can now be used to show easily that 
the codimension 1 duals always have positive net length. 
This is the content of the next result. 



Theorem 3 (Codimension 1) Let r be a codimension 1 
shared face of two n-dimensional simplices embedded in WL N , 
N > n forming a Delaunay pair. Then the signed length * t 
is positive. 

PROOF. When N — n, the results directly follows from 
Lemma 2 since the circumcenters are in the correct order. 
For N > n, we can isometrically embed the simplices in E" 
in which case, the circumcenters are again in the correct 
order and the result follows. 

In the N > n case, the signs of the elementary dual edges 
of * t are assigned in the affine spaces of the corresponding 
n-dimensional simplices. For example, consider a pair of 
triangles embedded in K 3 and meeting at a shared edge at 
an angle other than it. In this case, the signed length of the 
dual edge of the shared edge is determined as the sum of 
the two elementary dual edges which are measured in the 
planes of the two triangles individually. 

3.1. Dual of a Vertex in Triangle Mesh Surface 

Now we show that the area of the dual of an internal 
vertex in a pairwise Delaunay triangle mesh is always posi- 
tive. We prove this below by showing that the net dual area 
corresponding to a pair of triangles is positive. 
Theorem 4 Let r be an internal vertex in a pairwise De- 
launay triangle mesh embedded in M. N , N — 2,3. Then the 
signed area of-kr is a positive. 

PROOF, -kr is the Voronoi cell of vertex r in the pair- 
wise Delaunay mesh. Consider a pair of triangles sharing a 
common edge incident to r and if they are embedded in M 3 , 
isometrically project to M 2 (i.e., treat the shared edge as 
a hinge, and flatten the pair.) The circumcenters of these 
two triangles are in correct order by Lemma 2 and there 
are three possible cases as shown in Figure 4. Thus the net 
area of the two elementary dual simplices is positive when 
the signs are assigned using the rule described in Section 2. 
Summing over all edges containing r yields the full *r as 
a positive area. 



3.2. Dual of an Edge in Tetrahedral Mesh 

Theorem 5 Let t be an internal edge in a tetrahedral De- 
launay triangulation embedded in R 3 . Then -kr is a simple, 
planar, convex polygon whose signed area is positive. 

PROOF. *r of an internal edge r in a Delaunay triangu- 
lation may or may not intersect r. The vertices of *r are 
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Fig. 4. Elementary dual simplices of a vertex in a pair of triangles 
sharing an edge. The cases shown correspond to various positions of 
the circumcenters of the shared edge and the two triangles. 

circumcenters of tetrahedra incident to r and the bound- 
ary edges of *r are dual edges of triangles incident to r. 
Note that * r is the interface between the Voronoi cells cor- 
responding to the two vertices of r and thus is a bounding 
face of both Voronoi cells. Since the Voronoi cell of a vertex 
is a convex polyhedron [6], *r is simple, planar and convex. 

Suppose t intersects *r. Then the tetrahedra incident 
to r and the edges of *r have to be in a configuration 
shown in left part of Figure 5. A configuration in which the 
triangles incident to r are reflected about r is impossible 
due to Lemma 2. 

Now, to see that the signed area of * r is positive, consider 
two elementary dual simplices of * r incident to a shared 
face a of two tetrahedra in the fan of tetrahedra incident 
to t. These two elementary dual simplices can be in one of 
the two configurations as shown in Figure 6. In both cases, 
c T is the circumccntcr of the edge r, c CT is the circumcenter 
of the shared face cr, and c p and c\ are the circumcenters of 
the two tetrahedra. Also, in both cases, using the sign rule 
of Section 2 the sum of the signed areas of the elementary 
dual simplices is positive, and hence, the signed area of *r 
composed of these elementary dual simplices is positive. 

Next consider the case in which r does not intersect * r 
as shown in right part of Figure 5. A boundary edge of *r 
is called near side if it is visible from the midpoint of t, 
otherwise, it is called a far side edge. Figure 6 shows the 
net dual simplices of a near side and far side boundary edge 
of *t. By the sign rule of Section 2, far side elementary 



dual simplices have a net positive signed area while near 
side elementary dual simplices have a net negative signed 
area. The negative areas of the near side dual simplices are 
covered by the positive areas of the far side dual simplices. 
Thus, the sum of all these elementary dual simplices which 
is the signed area of *r is positive. 




Fig. 5. An internal edge r of a tetrahedral mesh may or may not 
intersect *t. The views here are along r which appears as a point. 
The short lines are half-planes of the triangles incident to r. The 
tetrahedra are labeled a, 6, etc. Each boundary edge of *r corre- 
sponds to the triangle indicated by the coloring. The half planes 
could potentially be a reflection about r but that is impossible in a 
Delaunay mesh due to Lemma 2. 




Fig. 6. Representative elementary dual simplices of * r when it inter- 
sects r (left side) and does not intersect r (right side) corresponding 
to the two cases shown in Figure 5. 

3.3. Dual of a Vertex in Tetrahedral Mesh 

Theorem 6 Let r be an internal vertex of a tetrahedra 
Delaunay mesh embedded in M 3 . Then the volume of-kr is 
positive. 

PROOF. * r of a vertex r in a Delaunay tetrahedral mesh 
is a convex polyhedron that is the Voronoi dual cell of r [6] 
and thus r is inside *r. The faces of *r are duals of edges 
incident to r. By Theorem 5 all these faces have a positive 
signed area. The direction corresponding to traversal from 
r to an edge center always has a positive sign. Thus each 



5 



pyramid formed by r and a boundary face of *t has positive 
volume. Thus, the volume of *r is positive. 

4. Requirements on Boundary Simplices 

In the previous section we have only considered internal 
simplices in a pairwise Delaunay mesh. For simplices lying 
in the boundary of a domain we require an assumption to 
ensure positive duals. We call a simplex a one-sided with 
respect to a codimension 1 face r if its circumcenter c a lies 
in the same half space as the apex with respect to r in the 
afhne space of a. 

We show below that the only assumption then needed is 
that a top dimensional simplex with a codimension 1 face 
in the domain boundary should be one-sided with respect 
to the boundary face. 

Consider a pairwise Delaunay mesh of dimension n em- 
bedded in M. N , N > n. Assume that r is an (n — 1)- 
dimensional face appearing in domain boundary and r -< 
a - ™ such that a 11 is one-sided with respect to r. 
Theorem 7 For a mesh such as above, a dual of codimen- 
sion 1 faces has positive length. For n — 2 and N = 2 or 3, 
and for n = N = 3, duals of all simplices at all dimensions 
have positive areas or volumes. 

PROOF. The codimension 1 dual of r in all cases has 
positive length using our sign rule since o~ n is one-sided with 
respect to t. As a result, for a surface triangle mesh, that 
is n — 2 and N — 2 or 3, it easily follows from our sign rule 
in Section 2 that the dual of a vertex on the boundary also 
has a positive area. 

For n = N = 3, one configuration for the dual of an edge 
r incident to the boundary is shown in Figure 7. In this fig- 
ure, the plane containing the codimension 1 faces incident 
to r are shown as short line segments, and the coloring of 
boundary edges of *r show the corresponding codimension 
1 face they are dual to. The other configuration in which 
the planes containing the faces incident to r are mirror im- 
ages of ones shown is not possible since then the circum- 
centers of tetrahedra will not be in the correct order as in 
Lemma 2. Thus, by our sign rule, all elementary dual sim- 
plices of -kr are positive and hence the signed area of *t is 
positive. Finally, it follows from our sign rule that the dual 
of a vertex on the boundary is also positive since each of 
the elementary dual pyramids will have a positive volume. 



5. Conclusions and Outlook 

For planar triangle meshes and for tetrahedral meshes in 
three dimensional space the condition of being pairwise De- 
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Fig. 7. Dual of an edge r lying in the boundary of a Delaunay 
tetrahedral mesh. The meaning of colors and small lines is as in 
Figure 5. 

launay is equivalent to being Delaunay. Thus most commer- 
cial and freely available meshing software can generate such 
meshes. In our experience, several codes for planar mesh- 
ing also generate meshes for which the one-sidedness condi- 
tion on the boundary is satisfied. For example, the popular 
meshing code called Triangle has an option for conforming 
Delaunay triangulation which generates Delaunay meshes 
with one-sided boundary triangles. For tetrahedral meshes 
with acute input angles this property may be harder to 
achieve. In general however, algorithms for creating tetra- 
hedral meshes with one-sided boundary tetrahedra do ex- 
ist [3,5,7]. Note that one-sidedness is equivalent to an "ori- 
ented" Gabriel property (using diametral half-balls) for the 
boundary faces. 

The pairwise Delaunay condition also appears to be more 
natural for DEC than other conditions that are used in 
place of Delaunay in the case of surfaces. For example, 
some researchers require that the equatorial balls of trian- 
gles not contain another vertex. This disqualifies surfaces 
with many folds or sharp turns. Another alternative is to 
define intrinsic Delaunay condition based on geodesies on 
the triangle mesh but algorithms for such surfaces can be 
complicated to implement. Yet another alternative is to use 
Hodge-optimized triangulations [9]. But creation of these 
requires an additional optimization step. On the other hand 
Hodge-optimized triangulation is a very interesting gener- 
alization of Voronoi-Delaunay duality with many applica- 
tions. 

The invention of algorithms that generate pairwise De- 
launay surface meshes is left for future work. So is the proof 
of our conjecture that the case of codimension other than 1 
has positive volume for general dimension and embedding 
space for pairwise Delaunay meshes with one-sided bound- 
ary simplices. Nevertheless, the practically important cases 
have all been settled by this paper. 
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Fig. A.l. See proof of Lemma 8 below. 



In fact here we prove a stronger result than Lemma 2. 
We will show that the circumcenters are in the correct or- 
der if and only if the pair of simplices is a non-degenerate 
Delaunay pair. 

Lemma 8 Let r be an (n — 1) -dimensional simplex in K™. 
Let L and R be points separated by r. Let c\ and c p be the 
circumcenters of the n-dimensional simplices A = L * r and 
p = R * t , respectively. Then, c\ and c p have the same 
relative ordering with respect to r as L and R if and only if 
A and p are a pair of non- degenerate Delaunay simplices. 



PROOF. Since A and p are a Delaunay pair, the afhne 
space of r separates L and R. See Figure A.l. Let c T and r T 
be the circumcenter and the circumradius of t, respectively. 
Now, c\ and c p lie on a line £ that passes through c T and 
is orthogonal to the affine space of r. Let h\ be the signed 
distance along £ from c\ to c T . Similarly, let h p be the signed 
distance from c p to c T . For now, it is sufficient that these 
distances be signed and whether the positive direction is 
along L or R is not important. Next, orthogonally project 
R onto £, and let r R be the (positive) distance from R to 
its projection onto £. Finally, let h R be the signed distance 
(along £) from the projection of R onto £ to c T . Notice that 
ha is necessarily either negative or positive depending on 
the choice of positive direction to be either along L or R, 
respectively. 

By elementary geometry, the squared circumradius of A 
is h\ + r 2 and the squared circumradius of p is h 2 p + r 2 . 
Similarly, the squared distance from c\ to R is (h\ — hji) 2 + 
r 2 R . Since A and p form a Delaunay pair, R lies outside the 
circumsphcre of A. Thus, the squared circumradius of A is 
less than the squared distance from c\ to R: 

h 2 x + r 2 T < (h\ - h R ) 2 + r% , 
=*> r 2 < r\ + h% - 2 h R h x . 

Also, since R lies on the circumsphere of p, the distance 
from c p to R is the same as the distance from c p to a vertex 
of r. Thus, we have: 

h 2 p + rl = {h p -h R f + r R , 

=> r r = r R + h 2 R -2 h R h p . 

Using this in the previous inequality, we obtain: 

r R + h R~ 2 h R hp <r 2 R + h 2 R -2h R h x , 
h R h p > h R h x ■ 

Finally, we choose a coordinate direction along £ to fix signs 
in the signed distances along £. If we choose the direction 
towards the half space containing R to be positive, h R is 
positive. (We will call this the positive i?-direction.) As a 
result, the last inequality above simplifies to h p > h\. This 
means that h p is larger along the positive indirection. If we 
choose the direction along L to be positive, h R is negative 
and we obtain h p < h\. In this case, h\ is larger along the 
positive L-direction. 

Conversely, if A and p are not a Delaunay pair, then the 
distance from c\ to R is less than the circumradius of A. 
Thus, all inequalities will reverse directions and therefore 
the circumradii will be in the wrong order. 
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